Emergent Photons and New Transitions in the 0(3) Sigma Model with Hedgehog 
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We study the effect of hedgehog suppression in the 0(3) sigma model in D=2+l. We show 
via Monte Carlo simulations that the sigma model can be disordered while effectively forbidding 
these point topological defects. The resulting paramagnetic state has gauge charged matter with 
half-integer spin (spinous) and also an emergent gauge field (photons), whose existence is explicitly 
demonstrated. Hence, this is an explicit realization of fractionalization in a model with global SU(2) 
symmetry. The zero temperature ordering transition from this phase is found to be continuous but 
distinct from the regular Heisenberg ordering transition. We propose that these phases and this 
phase transition are captured by the noncompact CP^ model, which contains a pair of bosonic 
fields coupled to a noncompact (7(1) gauge field. Direct simulation of the transition in this model 
yields critical exponents that support this claim. The easy-plane limit of this model also displays 
a continuous zero temperature ordering transition, which has the remarkable property of being 
self-dual. The presence of emergent gauge charge and hence Coulomb interactions is evidenced 
by the presence of a finite temperature Kosterlitz-Thouless transition associated with the thermal 
ionization of the gauge charged spinous. Generalization to higher dimensions and the effects of 
nonzero hedgehog fugacity are discussed. 



I. INTRODUCTION 

Since the initial proposal of fractionalization (phases 
where the elementary excitations are fractions of the 
electron) as the underlying explanation for the unusual 
properties of the cuprate superconductors jlj], there has 
been much progress in the theoretical understanding of 
these phases. Whether such fractionalized phases in the 
absence of magnetic fields in spatial dimension greater 
than one are actually realized in any experimental sys- 
tem is a matter of current debate. Fractionalized phases 
can exhibit properties that are strikingly different from 
conventional phases of matter, hence they are attractive 
candidates for modeling strongly correlated systems that 
exhibit anomalous behaviour. However, unambiguous ex- 
perimental evidence for the presence of such phases in any 
experimental system is still lacking. In part this may be 
because the correlations in such phases are subtle, and 
hence definitive experimental signatures are hard to de- 
vise. This provides a strong motivation to seek a deeper 
understanding of these phases. Furthermore, fractional- 
ized states have been proposed as a means to build quan- 
tum bits that are inherently robust against decoherence 
0. 

An important theoretical development has been the 
discovery of a number of microscopic models H, IE 
S II II H ES 113 that can be shown to exhibit this 
exotic physics. Although these microscopic models are 
defined in terms of bosons (or spins) on a lattice with 
short ranged interactions, fractional excitations that are 
charged under an emergent gauge field, as well as ex- 
citations of this emergent gauge field are obtained on 
solving these models. However till date, there have been 
no microscopic models available with the full SU(2) spin 
rotation symmetry. Indeed it is important to verify that 
the additional constraints imposed by spin rotation sym- 



metry do not exclude the possibility of fractionalization. 
Here we will describe a model that possesses the full spin 
rotation symmetry, but can be explicitly shown to exhibit 
fractionalization and possess an emergent gauge field in 
the deconfined phase. Moreover this model is found to 
have a quantum critical point with full spin rotation sym- 
metry, but which is distinct from the Heisenberg transi- 
tion. The properties of the fractionalized phase and the 
transition, as well as various deformations on the model, 
will be studied in detail in this paper. 

Most of the earlier work that constructed models ex- 
hibiting fractionalization engineered the energetics so as 
to select a low energy manifold. Constraining states to 
lie in this manifold introduces the gauge fields, which 
then need to be in the deconfined phase for fractionaliza- 
tion to occur. Here we will rely on a different route to 
fractionalization, which may be described as fractional- 
ization from defect suppression. Indeed this approach is 
closer in spirit to where the Z2 fractionalized state 
was regarded as a quantum disordered superconductor, 
obtained by proliferating even winding-number vortices 
while suppressing vortices of odd winding number. 

Here we will mainly be concerned with the 0(3) sigma 
model in D=2-fl. This model consists of 0(3) quan- 
tum rotors represented by unit three-vectors ('spins'), 
defined on the sites of a two dimensional spatial lattice. 
Neighbouring rotors are coupled via a ferromagnetic in- 
teraction. By the usual Quantum to Classical Statisti- 
cal Physics mapping, the ground state properties of this 
model can be conveniently mapped onto the physics of 
the Heisenberg model at finite temperature in three di- 
mensions. Clearly, there exist point topological defects 
in the three dimensional Heisenberg model that carry an 
integer topological charge, which simply correspond to 
hedgehog configurations of the spins. In terms of the 
quantum model, these are events in space-time (instan- 
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tons) which change the skyrmion number of the system. 
We now ask the question: Is it possible to disorder the 
three dimensional Heisenberg model in the effective ab- 
sence of the hedgehog defects? This has been a long 
standing issue, discussed in several works, for example 
[T^lIT^. but had not been conclusively settled. Here we 
will present fresh results from Monte Carlo simulations 
and arguments that convincingly demonstrate that the 
answer to this question is yes. 

Having established this, we will ask the question, 
what is the nature of this hedgehog-free paramagnetic 
phase P*? Although the spin-spin correlations are short 
ranged, the absence of hedgehog fluctuations lead to a 
'hidden order' in this phase. Indeed, it will be proposed 
that the physics of the hedgehog-free model is captured 
by what we will call the non-compact CP^(NCCP^) [l|| 
model in D=2-|-l. This model consists of a doublet of 
bosonic fields ('spinous') that transforms as a spinor un- 
der spin rotations, coupled to a non-compact U(l) gauge 
field ('photon'). In this representation, the hedgehogs 
correspond to the monopoles of the U(l) gauge field, 
and eliminating hedgehogs leads to the noncompactness 
of the gauge field. The NCCP^model has two obvious 
phases, one where the spinous are condensed which is 
the ferromagnetic phase, and the other where the spinous 
are gapped, which corresponds to an exotic paramag- 
net, with a gapless photon excitation. To verify that the 
paramagnetic phase P* obtained by the hedgehog-free 
disordering of the 0(3) sigma model is indeed the same 
phase as the paramagnet in the NCCP^model, we use the 
Monte Carlo method to measure the correlations of the 
spin chirality, which is roughly rii • (n2 x n^) for a triangu- 
lar face with three spins ni. 2, 3 • This should be equivalent 
to the flux correlations of the noncompact gauge theory. 
Indeed the expected long range correlation functions with 
the very characteristic dipolar form are found. Further- 
more, this correspondence implies that the ordering tran- 
sition of hedgehog- free 0(3) model is not in the Heisen- 
berg universality class, but in the same universality class 
as the ordering transition of the NCCP^model. This may 
be checked by comparing the universal critical exponents 
in the two models. Indeed, we find that the hedgehog-free 
0(3) model undergoes a continuous ordering transition 
with critical exponents that are clearly distinct from the 
Heisenberg exponents, but which are consistent with the 
exponents obtained in direct Monte Carlo simulation of 
the transition in the NCCP^model. This provides further 
evidence that the NCCP^ model captures the physics of 
the hedgehog-free 0(3) model. These exponents also turn 
out to be consistent with those obtained from an earlier 
attempt to disorder the hedgehog- free 0(3) model 

Given the central role played by the NCCP^model, 
we also present analytical results on some of the striking 
properties of this model. First, we consider the easy- 
plane deformation of the NCCP^ model when the full 
spin rotation symmetry is broken down to U(l) by the 
presence of easy-plane anisotropy. This model is found 
to have the amazing property that under the standard 



duality transformation, it maps back onto itself. In par- 
ticular the zero temperature ordering transition is found 
to be self-dual. Second, we consider the effect of finite 
temperature in the NCCP^model. We argue that there 
is a thermal Kosterlitz-Thouless phase transition out of 
the P* phase. This can be understood as an ionization 
transition of the logarithmically interacting spinons. The 
logarithmic interaction is of course the Coulomb poten- 
tial in two spatial dimensions and hence the presence 
of such a transition is proof of the existence of gauge 
charged particles. 

We now briefly comment on the relation of the present 
paper to earlier relevant work. Lau and Dasgupta (131] . 
considered the 0(3) model in three Euclidean dimensions 
on a cubic lattice and applied complete monopole sup- 
pression at every cube of the lattice. This strong con- 
straint led to the model always being in the ferromagnetic 
phase, and the exotic paramagnet was not uncovered in 
that work. Subsequently, in an important extension, Ka- 
mal and Murthy [ij allowed for a more flexible defini- 
tion of the no monopole constraint by allowing monopole- 
antimonopole pair fluctuations if they occurred on neigh- 
bouring cubes. In this way, they were able to obtain a 
disordered phase, and also found a continuous ordering 
transition with non-Heiscnberg exponents. However, as 
pointed out in that work itself, there are unsatisfactory 
features of this algorithm when closely spaced monopoles 
occur. The problem arises in the presence of loops of 
monopole and antimonopoles, where combining them in 
pairs is ambiguous. This allows for Monte Carlo moves 
that annihilate a monopole and an antimonopole belong- 
ing to different pairs, since the remaining monopoles and 
anti- monopoles on the loop can be 're-paired'. This re- 
quires making a time consuming nonlocal check for re- 
pairing, every time such an event is generated. Conse- 
quently, Kamal and Murthy had to resort to the approx- 
imation of making only local checks of re-pairing. The 
algorithm then also has the undesirable property that a 
given configuration of spins may be allowed or not al- 
lowed depending on the history of how it is generated. 
Here we will adopt a more restrictive condition - that 
monopole anti-monopole pairs are allowed only if they 
are isolated - to completely circumvent these problems. 
An improved procedure for defining the monopole num- 
ber allows us to easily work with more complicated lat- 
tices. This flexibility will prove very useful in obtaining 
the disordered phase. Furthermore, our identification of 
the NCCP^ model to describe this physics allows us to 
bring a whole series of tests to bear on the P* phase 
and the transition, which was not done previously in the 
absence of such an understanding. 

In this paper we will be completely suppressing the 
free hedgehog defects, i.e. setting their fugacity to zero. 
We now briefly discuss the effect of a finite fugacity, and 
possible relevance of the physics described here to sys- 
tems that may be realized in nature. It is well known 
that for pure U{1) gauge theories in D=2-|-l, introduc- 
tion of monopoles leads to confinement This result 
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also implies that the P* phase is unstable to the intro- 
duction of a nonzero hedgehog (or monopole) fugacity. 
However, if this fugacity is small to begin with, decon- 
fincment physics would play an important role in the fi- 
nite temperature or short wavelength properties of the 
system. 

While a finite monopole fugacity necessarily destroys 
the zero temperature P* phase, its effect on the phase 
transition, where critical gauge charged bosons are 
present which act to hinder monopole tunneling, is less 
obvious. In fact, in the sigma model description of the 
spin half antiferromagnet on the square lattice. Berry 
phase effects lead to a quadrupling of the monopoles 
[18], which are then more likely to be irrelevant as com- 
pared to single monopole insertions. If monopoles are 
then irrelevant at the critical point, the transition could 
still be controlled by the monopole suppressed 0(3) crit- 
ical point (or equally the critical NCCP^model). How- 
ever, monopole relevance in the adjoining P* phase im- 
plies that this transition is sandwiched between two con- 
ventional (not deconfined) phases. This dangerously- 
irrelevant-monopoles scenario is advocated in a forth- 
coming paper 19], which argues that it is possible to 
have a continuous transition between a Neel and Va- 
lence Bond Solid (VBS) state for the square lattice spin 
half antiferromagnet, which is in the same universality 
class as the hedgehog free 0(3) transition (i.e. the crit- 
ical NCCP^theory) studied here. In fact, Ref. ^'i'" 
gues that the easy plane version of this transition may 
already have been seen in the numerical experiments of 
|2C| . where there appears to be a continuous transition 
between a spin ordered state and a VBS in a square lat- 
tice spin half model with easy plane anisotropy. Such 
a direct transition is very natural in the dangerously- 
irrelevant-monopoles scenario, and would be controlled 
by the critical NCCP^ model with easy plane anisotropy, 
which is also studied in this paper. Thus, the transi- 
tions of the hedgehog suppressed 0(3) model, in both 
the isotropic and easy plane limits, could be realized in 
these situations even without explicit hedgehog suppres- 
sion, and might potentially be seen directly in nature. 
Finally, we note that in D=3+l a deconfined phase, with 
photons and gapped spin half particles, can in principle 
exist even with a nonzero hedgehog fugacity. 

The layout of this paper is as follows. In Section ^ 
we study the 0(3) sigma model with hedgehog suppres- 
sion. We begin by describing the particular lattice geom- 
etry and hedgehog suppression scheme used in the Monte 
Carlo calculations. We then present Monte Carlo results 
that show the presence of a spin disordered phase in our 
hedgehog suppressed model. We argue that the physics 
of the hedgehog suppressed sigma model is captured by 
the NCCP^model, which implies the presence of photons 
in this disordered phase. This leads to the prediction that 
spin chirality correlations in the disordered phase take on 
a very particular long ranged form which is tested in the 
Monte Carlo calculations. Next, we turn to the univer- 
sal properties of the ordering transition in this hedgehog 



suppressed model and find exponents that are distinct 
from the Heisenberg exponents. These are then com- 
pared with exponents calculated from directly simulating 
the NCCP^model. In Section IIIII we consider various 
deformations of the NCCP^ model that correspond to an 
easy plane anisotropy, a Zeeman field, and the effect of fi- 
nite temperature. In particular we prove the remarkable 
self duality of the easy plane model. Finally, in Section 
IIVI we briefiy discuss possible higher dimensional exten- 
sions of this physics. 

II. THE HEDGEHOG FREE 0(3) SIGMA 
MODEL 

A. The Model 

We perform Monte Carlo simulations of the three di- 
mensional classical 0(3) sigma model with hedgehog sup- 
pression. The lattice that we consider is a decorated cu- 
bic lattice as shown in Fig. with unit vectors iij at 
the vertices and edge-centers of the cubic lattice. As 
described below, monopole numbers are associated with 
the centers of each of these cubes. This choice allows 
for more spin fiuctuations when hedgehog suppression is 
applied than the simple cubic geometry of Refs. ITsl [3 
Neighbouring spins are coupled via the usual ferromag- 
netic Heisenberg interaction and hence allowed states are 
weighed with the factor e~^, with the energy function E 
given by: 

i? = -J^n,.n, (1) 

In order to define the monopole number in each cube, 
we follow reference pl| and first introduce an auxiliary 
variable, the gauge potential Aij between any pair of 
neighbouring sites with spin orientations n^, n^. This 
is defined by introducing an arbitrary reference vector 
and forming the spherical triangle (n*, n^, n^). The 
edges of a spherical triangle are of course segments of 
great circles. If the solid angle subtended by this spheri- 
cal triangle is 0[n*, n^, rij], then we define: 

1 + ■ rij + n» • + rij ■ rij + in* • (n, x n^) 
^J2{l -f n* • nj)(l -I- n* • nj){l + riirij) 

A different choice of the reference vector, n'^, only leads 
to a 'gauge' transformation of A: Aij Aij + Xi ~ Xji 
where Xi = n*, n^], etc. Thus, gauge invariant 

quantities are independent of the choice of the reference 
vector. Note also that e*-^"' — e~*-^^' and these gauge 
fields are only defined modulo 27r. Thus we have de- 
fined a compact gauge field in terms of the spins. We 
then define a flux Pg on every face bounded by the sites 
(l,2,...,n,l): 
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FIG. 1: a) Decorated cubic lattice used in the simulations. 
Spins live on the lattice points shown, and the monopole 
number is defined within each cube, b) The only spin config- 
urations accepted in the simulation are those that are either 
hedgehog free, or have hedgehogs that can be paired uniquely 
into isolated nearest neighbor hedgehog-antihedgehog pairs. 
A schematic depiction of such a pairing is shown here in a 
vertical section through an isolated pair. 



gether the problems in the work of 0|, where such an 
isolation of pairs was not demanded. 

To summarize, the statistical ensemble is defined as 
follows. For each spin configuration, we determine the 
hedgehog numbers associated with each cube of the lat- 
tice. If this sample clears the constraint of no free hedge- 
hogs, (we mentioned two versions of this constraint, full 
suppression constraint and the isolated neutral pairs con- 
straint), then this configuration is allowed in the ensem- 
ble and is weighted with a relative probability determined 
by the energy function JQ). 

We simulate this ensemble {221] using single spin 
Metropolis updates in the restricted configuration space. 
The data presented below is taken for 20,000-200,000 
Metropolis steps per spin. 



with F\j G (— 7r,7r]. Clearly, the fiux is gauge invariant 
and hence independent of the choice of reference vector 
ri:*. The physical meaning of the flux is most readily 
appreciated by considering a triangular face with spins 
Hi, 2, 3, where it is approximately the spin chirality: 

sinFA ni • (n2 X ns). (4) 

The hedgehog number k enclosed in a given volume 
is then the net flux out of this volume ^ — 2-Kk, 
which is guaranteed to be an integer from the previous 
definitions. Note that the hedgehog number is simply 
some function of the spins on a given cube. This defi- 
nition is identical to the traditionally used definition of 
hedgehog number for volumes that are bounded by tri- 
angular faces. For more complicated geometries (like the 
one employed in this work) however it is a much more 
natural and powerful definition, since it does not rely on 
an arbitrary triangulation of the faces and can be quickly 
computed. 

We now consider whether a disordered phase may be 
obtained while suppressing the hedgehog configurations. 
This will favour the ferromagnetic state which is clearly 
free of hedgehogs; indeed with full hedgehog suppression 
on the simple cubic lattice ,13i] an ordered phase was 
found even at zero spin coupling. In the decorated lat- 
tice shown in Fig. full hedgehog suppression in each 
cube seems to give rise to a disordered state for small val- 
ues of the spin coupling J. However, in order to open a 
larger window of disordered phase, and obtain more solid 
evidence of disorder in the system sizes available, we will 
allow for hedgehog-antihedgehog fluctuations on nearest 
neighbour (face sharing) cubes. In contrast to Ref.Q, we 
will only allow for configurations with isolated hedgehog- 
antihedgehog pairs; in other words if a cube contains a 
hedgehog of strength g, it must contain a nearest neigh- 
bour cube with a hedgehog of strength —q and no hedge- 
hogs in all other nearest neighbour cubes. This is shown 
schematically in Fig. ^ and gives an unambiguous pre- 
scription for combining the hedgehog and antihedgehogs 
into isolated, neutral pairs, and allows us to avoid alto- 



B. The Disordered Phase 

We now discuss the results of the Monte Carlo simula- 
tion with hedgehog suppression. First, in the absence of 
any hedgehog suppression, the system is found to have 
the usual Heisenberg ordering transition at Jc.Hois ~ 
1.7. Implementing hedgehog suppression that only al- 
lows neutral, isolated pairs of hedgehogs to occur, gives 
a smaller but still sizeable region < J < 0.7 over 
which the system remains magnetically disordered. This 
can be seen in Fig. |21 where the magnetization per spin 
TO is plotted for varying system sizes with linear di- 
mension L = 6,8, 12, 16 (the total number of spins is 
-^spin = 4L"^). The magnetization per spin is seen to ap- 
proach zero with increasing system size, for small enough 
values of J. A more convincing demonstration is made 
in the inset, where we plot the product of m and the 
square root of the total number of spins. For disordered 
spins, the average magnetization per spin is expected to 

— 1/2 

decrease as N . Indeed, as seen in the figure, this 
situation is realized at least for J < 0.5. 

One may nevertheless worry if there is some other spin 
order, such as antiferromagnetic or spiral order, that is 
not detected by the above zero-momentum magnetiza- 
tion. The most direct evidence against any magnetic or- 
der is obtained from the spin-spin correlation, which is 
found to be ferromagnetic throughout and rather short- 
ranged. For J = this is shown in Fig. 13 and the spin 
correlation indeed decays very quickly, with the correla- 
tion length of order one half lattice spacing. The above 
does not mean that the spins are completely uncorre- 
lated, rather that their correlation is more subtle as we 
will see below. 

This completes the evidence for the presence of a mag- 
netically disordered phase P* with suppressed hedge- 
hogs. We now investigate the nature of this paramagnetic 
phase. 
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FIG. 2: Magnetization per spin, m = {|M|)/iVspin, with M = 
^^iii, as a function of J for different system sizes (we show 
the data for both sweep directions). Inset shows the product 

1/2 

"^-^spini magnetically disordered phase, we expect the 

measured {|M|) ~ N^^^-i i^'^^ completely uncorrelated spins, 
the numerical coefficient is close to 1) 




FIG. 3: Spin-spin correlations for J = 0.0, measured for spins 
at the vertices of the cubic lattice separated by a distance 
r along the z direction. The system size is L = 16, so the 
measurements are done for r < 8. Note the logarithmic scale 
for the vertical axis; the lower cutoff is roughly the limit of 
what can be reliably measured in our Monte Carlo. 



C. Emergent Photons in P* 

We argue below that the P* phase of the hedgehog 
suppressed 0(3) model is distinct from the regular para- 
magnetic phase P in the model without such suppres- 
sion. The sharp distinction arises from the fact that P* 
contains a low energy photon excitation. While this is 
best understood by rewriting the 0(3) sigma model in 
the CP^ representation, we first provide a heuristic argu- 
ment for why such a low energy excitation may appear 
before passing to this more complete explanation. 

At any given time slice, the spin configuration (now of 



spins in a plane) can be given a skyrmion number. It 
is easily seen that hedgehog events change the skyrmion 
number by unity. Therefore, suppressing hedgehogs im- 
plies that the skyrmion number is a conserved quantity. 
Thus, if jo is the skyrmion density and ji,2 are skyrmion 
currents, they satisfy the conservation law 9^j^t = with 
/i = 0, 1, 2. This condition may be solved by writing 
j^l = ^^ivaduO-a, in which the skyrmion current is identi- 
fied with the flux of a U(l) gauge field. A natural dy- 
namics would then be given by a Lagrangian C — jfj,j^, 
which would give rise to a linearly dispersing photon. 

To gain further insight into the nature of the P* phase, 
we use the CP^representation of the 0(3) sigma model. 
It is well known that the pure 0(3) sigma model (no 
hedgehog suppression) can be rewritten in terms of a pair 
of complex bosonic fields z = (zi 2:2)'^ that is minimally 
coupled to a compact gauge field. The fields z transform 
as spinors under spin rotations, and have unit magnitude 
z^^z = 1. The spin vector is given by the bilinear n = 
z^^crz (where a are the Pauli matrices), and the flux of 
the gauge field corresponds to the skyrmion density of 
the original spin variables. Compactness of the gauge 
field implies the existence of monopoles, which act as 
sources or sinks of the gauge fiux. These are then to be 
identified with the hedgehogs which change the skyrmion 
number when they occur. Clearly, this CP^model has 
two phases, one where the z particles are 'condensed' 
which is the ferromagnetic phase (since the gauge neutral 
unit vectors n acquire an expectation value) , and another 
where the z particles are gapped. The gapped phase is 
essentially equivalent, at low energies, to a pure compact 
gauge theory which is known to be confining in D=2-|-l 
|16l |. This we associate with the regular paramagnetic 
phase. 

We now turn to a description of the phases with full 
hedgehog suppression within the CP^representation. In- 
deed, given the identification of the hedgehog defects 
with the monopoles of the CP^ theory, hedgehog sup- 
pression implies monopole suppression. This is most di- 
rectly implemented by passing to a noncompact gauge 
field which is free of monoples. The distinction be- 
tween a compact and a noncompact gauge theory cannot 
be overemphasized here, since it underlies all the new 
physics obtained in this work. The Euclidean action for 
the noncompact CP^model on a lattice is given by: 

S'nccpi = - Y (^^r+Ae""^ + c.c.) + y ^(A x a f 

(5) 

where the lattice curl is the sum of the gauge potentials 
around a plaquette. This NCCP^ model has two obvious 
phases - first, a phase where the z particles are 'con- 
densed' which is the ferromagnetic phase Second, a 
phase where the z particles are gapped - this we identify 
with the paramagnetic phase P*. However, the noncom- 
pact nature of the gauge field implies that there will be 
a low energy photon excitation in this phase. Indeed, 
within the NCCP^model, the asymptotic correlation of 
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the flux C)j,i, = (F^(r)f^(O)) in this phase is simply gov- 
erned by the free propagation of the photon which leads 
to the characteristic dipolar form: 



(6) 



In particular, for two points separated along the z direc- 
tion, we would expect 



Cyy{z) « 



(7) 



where B is some numerical coefficient. 

This prediction of the emergence of a photon in the 
P* phase may be readily checked by using the definition 
Eq. (|21 of the flux in terms of the spins on a face (the 
spin chirality), and studying its correlations 



C^^ir) = (sinP^(r)sinF,(0)) 



(8) 



in the hedgehog suppressed sigma model (as usual, F^^ 
is the flux through a face perpendicular to fi). The re- 
sults are shown in Fig. 0] which was taken deep in the P* 
phase with J = 0. This figure shows chirality-chirality 
correlations for points separated along the z direction and 
corresponds to the prediction in Eqn. Q . Indeed the ex- 
pected 1/r^ falloff is reproduced, as well as the sign of 
correlations and their approximate relative magnitude. 
[For the data in Fig. |3| we find Czz{z) « -1.7 Cyy{z); 
the slight discrepancy is most likely because the scaling 
regime is not quite reached for the separation in several 
lattice spacings.] Thus, although the spin-spin correla- 
tion function is short ranged in P* , the spin chirality 
correlations have long ranged power law forms. This is a 
result of the hidden internal order present in the system 
that arises from the suppression of hedgehog defects. It is 
this internal order that gives rise to the coherent photon 
excitation. (We also note that this topological order sur- 
vives in a small applied Zeeman field E —t E — h ^ ■ , 
which supports our claim for the gapped spinous, — see 
Sec. HiTbI below and Ref.El) 



D. The Transition 

We now study details of the ordering transition in the 
hedgehog suppressed 0(3) model. We use standard finite 
size scaling analysis in order to estimate the correspond- 
ing critical indices. First, to find the critical point with 
good accuracy, we measure the cumulant ratio (Binder 
ratio) 



9 = 



(9) 



It is expected to have the finite size scaling form 

giJ,L)=g{SL'/n , (10) 



.|Q-i Chirality-chirality correlations 
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FIG. 4: Chirality-chirality (flux) correlations measured along 
the z direction for the same system as in Fig. |21 Note the 
logarithmic scale for both axes. We also show a ~ 1/z^ hne 
to indicate the observed power law fallofi^. 
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Scaling plots with v=1 .0, 3/v=0-80 
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FIG. 5: Finite-size scaling plots for the cumulant ratio (left 
vertical axis) and magnetization (right axis), corresponding to 
the scahng forms lHH and llT^ . We used Jc = 0.725, u = 1.0, 
and P/u = 0.80; the range of the horizontal axis corresponds 
roughly to J G [0.40, 1.05] for L = 8 (compare with Fig.|5J. 



where S = J — Jc is the deviation from the critical point. 
In particular, the curves g(J, L) for different fixed L all 
cross near the critical Jc, which we estimate to be Jc = 
0.725 ± 0.025. Using the above scaling form, we also 
estimate the correlation length exponent 

V = 1.0 ± 0.2 . (hedgehog suppressed 03) (11) 

The corresponding scaling plot is shown in Fig. [S] 

To estimate the exponent /3, we study the finite size 
scaling of the magnetization per spin 



m{J,L)^ L-f^'^fiSL^'") 



(12) 



Our best scaling of the data is also shown in Fig. |31 and 
we find 

(3/v = 0.80 ± 0.05 . (hedgehog suppressed 03) (13) 
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jM CCP scaling plots with v=0.80. p/v=0. 




(J-Jc) L 



FIG. 6: Finite-size scaling study of the ordering transition in 
the NCCPVodel Eq. 0. The system is at K = 0.6 and the 
transition from the P* to the ferromagnetically ordered phase 
is observed at Jc = 1.255±0.02; the measured order parameter 
is defined from the field = zlcrzr- The scaling analysis is 
similar to Fig.|S] The exhibited plots are for i/ = 0.80, f3/v — 
0.85; the horizontal range corresponds to J £ [1.10, 1.40] for 
the L = 8 system. 



We remark that these exponents are consistent with 
the exponents obtained in the study by Kamal and 
Murthy 0| on the same model but using a different 
monopole suppression scheme. These exponents are 
clearly different from the three dimensional Heisenberg 
exponents which are accurately known to be t'Hcis = 
0.705 ± 0.003 and (^/i^)Heis = 0.517 ± 0.002. ^ 

The specific heat exponent a = 2 ~ di^ is expected 
to be negative for this transition (d is the space-time 
dimensionality), a ~ —1.0, hence a cusp singularity is 
expected here. Although such singularities are harder 
to detect than a divergence, we nevertheless look for it 
in the Monte Carlo study of the hedgehog suppressed 
model. It is found however that the specific heat remains 
completely featureless across the transition; the reason 
for this behavior is unclear. 

An important check on the above results in the hedge- 
hog free 0(3) sigma model is provided by comparing with 
direct simulations on the noncompact CP^ model given 
by the Euclidean action (|SJ). A complete numerical phase 
diagram of this model is reported elsewhere psl ] , here we 
focus on the transition between the P* and F phases and 
comparing the exponents with those obtained from the 
hedgehog free 0(3) sigma model. Indeed we find that 
this transition is second order with exponents 

z^ = 0.8±0.1, (3/iy = OM±Om , (NCCP^) (14) 



where the exponent /3 of course describes the onset of 
ordering in the gauge neutral field n — z^crz. The cor- 
responding scaling plots for this transition are shown 
in Figure El Clearly, these are consistent with the ex- 
ponents we obtained earlier for the hedgehog free 0(3) 
sigma model, which leads us to conclude that this tran- 



sition is indeed distinct from the Heisenberg ordering 
transition and is instead described by the deconfined 
to Higgs transition in the NCCP-'^model. The remain- 
ing difference between our best estimates for the criti- 
cal indices [and also between the apparent values of the 
Binder cumulant at the critical points g{Jc) ~ 1.35 for 
the NCCP-'^transition, and g{Jc) ~ 1.25 for the hedge- 
hog suppressed 0(3) transition] we attribute to the small 
system sizes considered. For example, for the sizes stud- 
ied here, the raw data (not shown here) Binder cumulant 
crossing in the NCCP^ model has a clear downward trend 
with increasing system sizes, while it changes more slowly 
in the hedgehog suppressed 0(3) model; this would nar- 
row the difference on going to larger system sizes. 

Finally let us note an important physical distinction 
between this transition and the Heisenberg transition, 
which is brought out by comparing the rj exponents. 
This exponent is related to the anomalous dimension of 
the n field and is given by = 2 — d + 2{/3/iy). For 
the Heisenberg transition, this exponent is very small 
1] — 0.033 ±0.004. However, for the hedgehog suppressed 
0(3) transition we find this exponent to be fairly large 
r] 0.6 (and rj « 0.7 from direct simulations of the 
NCCP^model). A large rj is to be expected if the magnon 
field n can decay into two deconfined spinous. Indeed, in 
the limit of noninteracting spinous, this exponent is ex- 
pected to approach unity. Hence, the large 77 obtained at 
this transition is to be expected on physical grounds. 



III. DEFORMATIONS OF THE NCCP^MODEL. 
EASY PLANE, ZEEMAN FIELD AND EFFECTS 
OF FINITE TEMPERATURE. 

In this section we consider various deformations of 
the NCCP^model, and the effect of finite temperatures. 
The motivation for this is twofold. First, understand- 
ing the behaviour of this model in these limits will pro- 
vide us with a whole slew of potential checks to further 
strengthen the identification between the NCCP^ model 
and the hedgehog suppressed 0(3) model. Some of these 
explicit checks, like the correspondence in the easy plane 
limits and the effect of a Zeeman field on the isotropic 
models, will be presented elsewhere [25j, while others are 
left for future work. Second, we will see below that these 
models have several interesting properties and may them- 
selves be directly realizable. For example, the continuous 
transition seen in the numerical experiments of |2C| has 
been conjectured in to be described by the easy plane 
NCCPi model. 



A. Easy Plane Deformation 

We first consider modeling the easy plane deformation 
of our 0(3) invariant sigma model with hedgehog sup- 
pression. This can be accomplished, e.g., by having fer- 
romagnetic interaction between neighboring spins of the 
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form — J(nfn^ + n^rij)- Clearly, in this case the spins 
would prefer to lie in the x — y plane in spin space; the 
global spin symmetry of this model is now broken down 
from 0(3) to 0(2) x Z2, where the 0(2) corresponds to 
spin rotations about the z axis in spin space, and Z2 
arises from symmetry under — > — n^. The appropri- 
ate deformation of the NCCP^ model involves applying 
the term J7(|2:i|^ -I- \z2\'^) at each site with [/ > that 
favours equal amplitudes for the two components of the 
spinor field. The corresponding n vector will then like 
to lie in the easy plane. We will call this the easy plane 
NCCPimodel. 

In fact it will be useful to consider the limit of extreme 
easy plane anisotropy, where the amplitude fluctuations 
of the z fields are frozen out, and only the phase fluctu- 
ations remain. The z field can then be parametrized by 

Zi(2) — ^ , and the direction of the spin in the easy 



plane is given by ~ 2zJ'z2 — e^^'*'^^'*''^' . The par- 

tition function on a lattice in three Euclidean dimensions 
can then be written as: 



Zep = 

Sep = 



/TT 
[i?0ii^(/.2]e"^^^ , (15) 
-TT 

-J ^ [cos(A^(/)i - a^) 4- cos(A^(^2 - a^)] 



r.fi 



K 



where the integration over the gauge fields is performed 
after suitable gauge fixing to ensure finiteness of the par- 
tition function (hence the prime in the integration mea- 
sure: [I?a]'). Below, we will consider applying the stan- 
dard duality transformations on this model [20. l27l | . For 
this purpose it is more convenient to pass to the Vil- 
lain representation which makes use of the approximate 

^2 

rewriting e"^^"^" « J2n'^^"°'^~'^ ■ Note that the Vil- 
lain form has the same 27r periodicity, a — > a -I- 27r, as 
the original partition weight, and for large J we will have 
J J . The Villain form is expected to retain the univer- 
sal properties and phase structure of the original model 
and therefore we will start with it and consider a series 
of exact transformations that perform the duality. The 
partition function written in the dual variables will be 
seen to be essentially the same as that written in the di- 
rect variables thus establishing the 'self-dual' nature of 
these models. 

The Villain form of the partition function is: 



-00 J —71 ■ • 



(Axa) 



J1,J2 



X e 2j 



VE(j?+j2)p*i:[(^</'i-a)ji + (A02-a)j 



^1, (16) 



where ji_ 2 are integer current fields that live on the links 
of the lattice. We first write Zy as a model in terms of 
the current loops only, and then perform an exact rewrit- 
ing (duality) in terms of the vortex loops. The two forms 



will have essentially identical characters when viewed as 
integer loop models; with proper identification of the cou- 
pling constants, this is the statement of self-duality. 

We begin by integrating out the phase fields 4>i,2- This 
imposes the condition that the integer current fields ji^2 
are divergence free. These integer loops are simply the 
world lines of the two bosons Zi_2, which carry the same 
gauge charge and hence their sum couples to the gauge 
field. Integrating out the gauge field gives rise to a long 
range Biot-Savart type interaction between the total cur- 
rents. The result can be written as: 



iih 

Xe^^(ji+j2)'-5+(i-,r')(ji+j2)p' 
Xe^^(ji^2)re-(r,r')(jl-j2),./ 

1 1 



Q+{k) 
G-ik) 



2 J 4K{sin^ !f + sin^ ^ 
1 

2J ■ 



sin' % 



(17) 

-(18) 

I 

(19) 



Thus, the original model is equivalent to a model of inte- 
ger current loops |28|, where the combinations ji—j 2 have 
short range interactions, while the gauge charged current 
combinations ji -|- j2 have long ranged interactions with 
the asymptotic form G+{R 00) = j^^- 

We now perform an exact duality transformation on 
((T7|l . First, the divergence free current loops can be 
written as the lattice curl of a vector field ji_2 = A x 
Ai_ 2/277 that lives on the links of the dual lattice, and 
the integer constraint is implemented using the identity: 



(2^rE„^(A-2™) = El 



ilA 



Integrating out the 



fields Ai^ 2 we obtain the following form for the partition 
function in terms of the dual integer current loops I1.2: 



Zv oc ^ (S(A • li)5(A • I2) 
I1.I2 



ll:l2 

^g-i(li+l2)..F+(r,r')(: 



xe 



:ii+i2),' 

■)(ll-l2),' 



J + 2if (sin' \ + sin' ^ 

sin'^) 



sin'^ 



^-^^'^ 2(sin'^ + sin'^ 



(20) 
-(21) 

I 

(22) 



One can argue that the two integer loops of this dual 
representation H2()(l are precisely interpreted as the vor- 
ticities in the original two boson fields. A vortex of equal 
strength in both boson fields (li = I2) can be screened 
by the gauge field, and hence has only short range in- 
teractions, while unbalanced vortices which cannot be 
completely screened interact via long range interactions. 
Correspondingly, the combination of the dual currents 
li -|- I2 interact via a short range interaction J-^{R), 
while the other combination li — I2 has a long range 
interaction with the same asymptotic form as in H18|l . 
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T- [R — >■ 06) = ^ . Indeed the partition function writ- 
ten in terms of dual loops (|20|l is essentially the same 
as when written in terms of the direct variables p7|l . 
which can be seen by making the association li <-> ji 
and I2 — j2- The only differences are in the form of 
the short range interactions which do not affect univer- 
sal properties. As argued below this immediately im- 
plies that the transition between P* and F phases in this 
model will be self-dual. Of course, we can also rewrite 
(Pn|l in terms of the fields conjugate to the currents li^ 2, 
obtaining the action in terms of the dual (vortex) fields 
and the dual gauge field, which will have essentially the 
same form as the original action Eq. (|16|l . We do not spell 
this out here since the exhibited forms already suffice. 

We first note the description of the various phases in 
terms of the properties of boson current loops and also 
the vortex (dual) loops. In the direct representation, 
when the partition function is dominated by small loops 
of ji ± j2, the system is in the 'insulating' or paramag- 
netic state P* . In this phase there is a single low energy 
excitation - the photon. On the other hand, when these 
loops become arbitrarily large, which certainly occurs if 
both J', if ;» 1, we are in the 'superfluid' or ferromag- 
netic phase F . Here too we have a single low lying mode, 
the magnon, which is the Goldstone mode arising from 
the spontaneously broken spin symmetry within the easy 
plane. 

The two integer loops of the dual representation H20() 
correspond to the world lines of the vortices in the two 
boson fields. It is easily seen, for example by analyzing 
the dual action in different parameter extremes, that the 
P* phase occurs if large vortex loops of both li ± I2 pro- 
liferate, while if both the vortex loops are typically small, 
the F phase results. 

In terms of the direct boson variables, the F phase 
is the ordered phase, and hence has a low energy Gold- 
stone mode, while the P* phase is the disordered phase 
and has a low energy photon mode. In the dual variables 
the roles are reversed - the F phase is the disordered 
phase with the dual photon, while the P* phase is the 
dual 'ordered' phase, with the Goldstone mode. In par- 
ticular, a direct transition between P* and F can either 
be thought of as an ordering transition in the direct vari- 
ables, or the reverse in terms of the dual variables. This 
interchange of ordered and disordered phases is typical 
of a duality transformation. What is special in this case 
is that the partition function is essentially identical when 
written in terms of the direct and dual variables. This 
implies that the transition must be self-dual. One direct 
consequence of this self duality is that the critical expo- 
nent (/3) measuring the rise of the order parameter on 
the ordered side of the transition is equal to that of a 
suitably defined disorder parameter on the opposite side 
of the transition. Further consequences arising from the 
self duality and other phases contained in this easy plane 
NCCP^model will be discussed at length elsewhere psj . 
along with results of Monte Carlo studies. Here, we just 
note that a continuous transition from the P* to F phase 




J 



FIG. 7: Effect of a Zeeman field on the zero temperature 
phases of the NCCP^model. (a) Isotropic case - Two distinct 
paramagnetic phases exist in the presence of a Zeeman field. 
The P* phase has a low lying excitation, the photon, while 
the P phase is gapped and breaks no global symmetries. The 
two are separated by a continuous transition in the 3D XY 
universality class (the arrow indicates the ordered to disor- 
dered 3D XY transition). At zero field the transition Ci is 
described by the critical NCCP^model. (b) Easy Plane case - 
Zeeman field is applied perpendicular to the easy plane. Spins 
can order in the easy plane to give the F phase which has a 
low lying excitation, the Goldstone mode. This is separated 
from the disordered P phase by a 3D XY transition. There 
are two paramagnetic phases, P" with a photon and a gapped 
P phase separated by a 3D XY transition. At zero Zeeman 
field the two 3D XY transitions meet to give the self dual 
CxY transition. 

is obtained, with critical indices: 

1/ = 0.60±0.05 , = 0.70±0.05 (easy-plane). (23) 

Finally, we note that this critical point is conjectured 
in [igj to control the continuous transition seen in the 
numerical experiments of . 

B. Effect of Zeeman Field 

We now consider applying a uniform Zeeman field 
along the direction in spin space, Sz — —h^^nl, 
and study the effect on the various phases and phase 
transitions in both the isotropic and easy plane models. 
In the CP ^representation this term can be written as: 

5z = -/^^(klp-|z2n (24) 

which breaks the zi ^ zi and the only remaining sym- 
metry is that of U(X) spin rotations about the axis. 

For the model with the isotropic exchange coupling, 
adding the term (|24|l to the action (O will result in 
the phase diagram shown in Figure |7^). At zero 
field we have the P* and F phases separated by the 
monopole suppressed 0(3) transition, or equivalently the 
NCCP^ordering transition C/. Turning up the Zeeman 
field, the ordered moment of F locks into the field di- 
rection - the resulting phase does not break any of the 
global symmetries and has no low lying modes. Hence 
we will call this phase P. (Note that for the isotropic 
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NCCP-'^model at zero temperature, this phase is only 
present when a Zeeman field is applied) . This phase can 
be understood as a Higgs phase with condensation of the 
zi field but not of the Z2 field. 

On the other hand, in the P* phase there is a gap to zi 
and Z2 particles. Turning up the Zeeman coupling lowers 
the zi branch until at the transition it condenses to give 
the P phase. Thus, the transition between the P* and P 
phases involves condensing a single scalar field coupled to 
a gauge field and hence is expected to be in the inverted- 
XY universality class (i.e. the P* to P transition is 
the ordered phase to disordered phase transition of the 
3D XY model). 

The shape of the phase boundary near the zero field 
transition can be related to the critical exponents of the 
NCCP^model. In particular it depends on the ratio of 
the scaling dimension of the Zeeman operator (y/j — d—^ 
for the isotropic case) to the scaling dimension of the 
'thermal' operator {yt = ■^). Thus, if SJ is the deviation 
away from the zero field critical point C/ , then the phase 
boundary is given by the curve 



plane model is essentially self dual, and this remains true 
on adding the Zeeman term. This explains the 'reflec- 
tion' symmetry of the universal properties in the phase 
diagram - for example the duality interchanges the P* 
and F phases while the P phase is taken to itself. This 
also implies that the transition from P* to P must be in 
the same universality class as the transition from F to 
P, which was indeed the result of the previous analysis, 
which finds them both to be in the 3D XY universality 
class. 

The shape of the phase boundaries near the zero field 
easy plane transition can be related to the critical ex- 
ponents of the easy plane NCCP^ model using Equation 
ESI While the thermal eigenvalue in this case is already 
known (Eqn. I23|l . the scaling dimension of the Zeeman 
operator needs to be determined. This is conveniently 
done within the easy plane NCCP^ model at criticality, 
by studying the scaling dimension of the operator that 
gives rise to unequal hopping strengths Ji =^ J2 for the 
two species of bosons in Eqn. 1161 To this end, consider 
defining the link operators: 



h, cx ((JJ)"^ (25) 
= Vh/yt (26) 

. The expressions for the scaling dimensions yield = 
Vh/yt ^ h'd — p. Using the numerical values of the expo- 
nents in Eq. (|14|) . we have <j> « 1.7, which gives a phase 
boundary that approaches the isotropic transition point 
horizontally as shown in Figure [7^. 

For the model with easy plane exchange coupling, the 
effect of adding the Zeeman field perpendicular to the 
easy plane can be argued to lead to the phase diagram in 
Figure \7}p. Here we avoid using the extreme easy plane 
limit, where was set identically to zero, in order to 
have a finite coupling to the applied field. (Alternatively, 
we can model the effect of the Zeeman field by considering 
an action similar to Eq. 1161 but with different couplings 
Ji and J2 for the two angle variables.) At zero field we 
have the P* and F phases separated by the monopole 
suppressed easy plane transition. Again, in the presence 
of a Zeeman field, a P phase is possible, with no sponta- 
neously broken symmetry and no low lying excitations, 
where only the zi field is condensed. The phase transition 
between P* and P is then driven by the condensation of 
the gauge charged zi field, and is hence expected to be in 
the same universality class as the ordered to disordered 
transition of the 3D XY model "2^. Starting in the P 
phase, as the exchange coupling is further increased, it 
becomes favourable for the Z2 field to also condense, so 
that there develops a finite expectation value for the in- 
plane spin operators + iny = 2z^Z2- This leads to 
the F phase, which spontaneously breaks the remaining 
spin rotation symmetry and has one Goldstone mode. 
Since the gauge field has already been 'Higgsed' away, 
this transition is the regular ordering transition of the 
3D XY model. As described previously, the original easy 



Ofjir) = {cos(A,^i - a) ± cos(A02 - a)} (27) 

Clearly, adding ^Ij '^i^l rise to unequal hop- 

pings, and hence has an overlap with the Zeeman oper- 
ator. Similarly, ^tj ^^^^ have an overlap with the 
'thermal' operator. Thus, studying the power law decay 
of the correlators of Of^ at criticality allows us to extract 
the scaling dimensions of the Zeeman and thermal oper- 
ators. The latter may be compared against other, more 
accurate determinations of the same quantity and serves 
as a check of this approach. Thus we find: 

yh = 1.2 ±0.3 (28) 
yt = 1.6 ±0.2 (29) 

the value of yt obtained is consistent with the value of v 
quoted in Eqn. 1231 from which yt ~ 1.7. Thus we have 
4> = 0.75 ± 0.2. Note that since < < 1, the shapes of 
the phase boundaries are as shown in Figure Ul^ |29|. 

The NCCP^ representation of the hedgehog- free easy- 
plane model in the presence of a Zeeman field also arises 
in the context of two component superconductors, as de- 
scribed in Ref. [30| where the continuum version of this 
identification was also noted. 

Finally we note that these properties of the easy plane 
model in the presence of a perpendicular Zeeman field 
arc potentially useful to testing the conjecture of 
that the continuous transition seen in the numerical ex- 
periments of [23 is controlled by the critical easy plane 
NCCP^model. An important point is that given the 
mapping of the Neel field to the unit vector of the sigma 
model, the Zeeman field considered here actually corre- 
sponds to a staggered Zeeman field on the spins of Ref. l20l 
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ISOTROPIC 




EASY PLANE 




(or noncompact CP") model will simplify the discussion 
and the very same results hold for the models of interest, 
i.e. the isotropic and easy-plane NCCP^ models, because 
it is only the logarithmic binding/unbinding of the gauge 
charged objects that matters. Since the finite temper- 
ature physics described above occurs in the charge and 
electric field sector, it will be useful to consider writing 
the thermal partition function Zji = Tr e~^^ in terms of 
the zero Matsubara frequency components of these fields. 
The appropriate effective energy is: 



T 

FIG. 8: Effect of finite temperature T on the NCCP Vodel. 
(a) Isotropic case: The thermal Coulomb phase Py, a finite 
temperature analogue of P* , has power law electric field corre- 
lations (Eq. 1331 and is separated from the more conventional 
P phase by a Kosterlitz-Thouless phase transition. This can 
be understood as the thermal ionization of spinons interacting 
via the 2D Coulomb (log) potential. The F phase is disor- 
dered at any finite temperature, (b) Easy Plane case: The 
ordered phase now has algebraic correlations at finite temper- 
ature (XY) which is also separated from the P phase by a 
KT transition. 



C. Finite Temperatures 

We now investigate the finite temperature properties 
of the NCCP^ model (or equivalently, the monopole sup- 
pressed 0(3) sigma model), both for the isotropic as well 
as the easy plane case. The main result is that there is 
a finite temperature version of the P* phase, which we 
call the thermal Coulomb phase P^, with power law cor- 
relations of the electric fields. This is distinct from the 
regular paramagnetic phase P with short ranged corre- 
lations, and is separated from it by a finite temperature 
Kosterlitz-Thouless phase transition. The existence of 
such a transition can be seen from the following physi- 
cal argument. In the zero temperature P* phase, there 
are gapped spinons that interact via a logarithmic inter- 
action. If this interaction persists to finite temperatures, 
then clearly for small enough temperatures, thermal fluc- 
tuations will only generate neutral spinon pairs. This is 
the P^ phase. However, at some higher temperature it 
becomes entropically favourable to proliferate unbound 
spinons (exactly as with logarithmically interacting vor- 
tices in the two dimensional XY model) and a Kosterlitz- 
Thouless type transition to a spinon-plasma phase will 
be expected. Indeed, for the hedgehog suppressed 0(3) 
model, the existence of such a finite temperature transi- 
tion may be viewed as evidence for the existence of emer- 
gent gauge charged objects interacting via a Coulomb 
potential, that happens to be logarithmic in two spatial 
dimensions which gives rise to the transition. 

We now sketch a derivation of these results - we choose 
to do this in a model with a single species of 'spinon' or 
gauge charged particle. Treating this single component 



i?{«r,e.4 = ^^n2 + ^^e2^ (30) 

r rfi 

where r and /i run over the sites and directions of the two 
dimensional spatial lattice. The partition function in- 
volves a sum over all configurations of the integer charge 
fields fir and real electric fields Cr^ (since we are working 
with a noncompact gauge theory) that satisfy the Gauss 
law constraint A • e = at each site of the lattice. Note 
that the boson number and electric field variables that we 
use here are conjugate to the phase and gauge potential 
variables that have been used so far (cf. Eas. I5ll6|l . Our 
treatment here is precise in the regime where the total 
Hamiltonian is dominated by the exhibited E{nr, Crfj,} 
part, but is also valid more generally in the effective 
sense. 

The Gauss law constraint can be implemented using a 
Lagrange multiplier ao to give the following expression 
for the partition function: 

{rir} 

(31) 

To establish that there is a thermal transition, we inte- 
grate out the fields ao(r), e^^, so the partition weight is 
now just a function of the integer charges. This yields: 

^e-^S^^'"-^"-'"-' (32) 

where for large separation the potential between charges 
takes the form ¥„' ~ j^k! -^'-'S \ r-r' \ ' which is the 
Coulomb interaction in two dimensions. It is well known 
that such a Coulomb gas has two phases, one with the 
charges bound into neutral pairs at low temperatures, 
and a plasma phase at high temperatures separated by a 
KT transition which occurs at Tc < {47tIC)~^. This is not 
surprising since we know that NCCP" model is dual to 
the 0(2) quantum rotor model in two spatial dimensions, 
and the spinons of the NCCP° map onto the vortices of 
the 0(2) rotor model, which interact logarithmically. 

A sharp distinction between the low and high tempera- 
ture paramagnetic phases can be drawn by looking at the 
electric field correlators. In the low (but nonzero) tem- 
perature P/^ phase, these fall off as the inverse square 
of the distance, while in the high temperature 'spinon 
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plasma' phase P, these correlators are short ranged: 

(e.Me.(0)) - {^iXa^gtr^ rtrl (33) 

where the indices fi, ly E {x,y}. As shown elsewhere 
|25l |. exactly the same results are obtained for the 
isotropic NCCP-'^model, and the phase diagram is shown 
in Fig.|2t- In particular, despite the presence of a global 
SU(2) symmetry, the transition remains KT. 

For the easy plane NCCP^ model, the expected finite 
temperature phase diagram is shown in Fig.lHb. Here, be- 
sides the discussed and P phases, there is also a power 
law phase (XY) that appears out of the zero-temperature 
ferromagnetically ordered phase and terminates at the 
usual KT transition to the P phase. Again, the similar- 
ity between the P^ and XY sides of this phase diagram 
have to do with the self duality of the underlying easy- 
plane NCCPi model. 



IV. TOWARDS D=3-|-l AND LAYERED 
PHASES 

All of the preceding discussion was focused on the 
D=2-t-l dimensional system where complete hedgehog 
suppression was required in order to obtain the decon- 
fined phase. It is natural to ask whether similar physics 
can be obtained in models with only a finite energy cost 
for hedgehogs. This can occur in the presence of criti- 
cal bosonic fields as argued in 19] leading to deconfined 
quantum criticality in D=2+l dimensional systems. An- 
other way in which a deconfined phase can be stabilized 
with a finite hedgehog fugacity is to consider D=3-|-l di- 
mensional systems. Now, the hedgehog is a particle, and 
a finite hedgehog fugacity does not necessarily destabi- 
lize the Coulomb phase. Thus, in principle one could ob- 
tain a D=3-|-l U(l) deconfined phase with global SU(2) 
symmetry, by disordering the 0(3) sigma model with a 
sufficiently large but finite hedgehog core energy. This 
question is currently under investigation; the obvious ex- 
tensions of the presented D=2-|-l dimensional realization 
either led to ferromagnetic order or to the conventional 
paramagnet. 

An interesting possibility that we remark upon is a lay- 
ered system that lives in three spatial dimensions, where 
each layer realizes the D=2-|-l dimensional deconfined 
phase. Such a phase exhibits quantum confinement in 
that the gauge charged spinous can move freely within 
a layer but not between the layers, and the photon re- 
mains a D=2-|-l dimensional particle. This is reminiscent 
of the sliding Luttinger liquid phase [sj of coupled one 
dimensional systems that also exhibits the phenomenon 
of quantum confinement, and has a low energy phonon 
excitation (which is the analogue of the photon in this 
lower dimensional system). 



V. CONCLUSIONS 

We conclude by highlighting the main results of this 
work. We considered the question of whether the 0(3) 
sigma model can be disordered without monopoles, and 
answered it in affirmative. We provided an explicit ex- 
ample of an 0(3) Heisenberg spin system with no free 
hedgehogs and no magnetic order in (2-1- 1)D. While 
this is in agreement with the earlier work of Kamal 
and Murthy^^, the potentially problematic features of 
the procedure adopted in that work have been entirely 
avoided in the present paper. 

Furthermore, we identified the proper description of 
the hedgehog suppressed 0(3) model which involves 
spinous coupled to a noncompact gauge field, which may 
be called the noncompact CP^model. Thus the hidden 
topological order resulting from hedgehog suppression 
gives rise to a photon-like low energy excitation in the 
paramagnetic (P*) phase of this model, which leads to 
power law correlations of the spin chirality. This may 
also be viewed as an example of a U(l) fractionalized 
phase (albeit with complete monopole suppression) with 
full SU(2) spin rotation symmetry. 

Our understanding of the hedgehog suppressed disor- 
dered sigma model is interesting from the statistical me- 
chanics point of view and addresses some long standing 
questions. In a sense, we identified how to "decompose" 
the 0(3) model into a part that involves the topologi- 
cal defects (hedgehogs), and the part which does not in- 
volve these. Such decomposition of the 0(2) model into 
the vortex and spin wave parts is well known. The cor- 
responding "spin wave" part for the 0(3) model turns 
out to be the NCCP^model, which clearly has little to 
do with spin waves. In particular, the commonly asked 
question whether the spin waves (perhaps non-linearly 
coupled) can disorder the 0(3) model seems to have no 
meaning. 

Another question has been on the role played by the 
hedgehog defects at the Heisenberg transition. A sharp 
formulation of this question is whether the ordering tran- 
sition in the hedgehog suppressed model is identical to 
the Heisenberg transition. Our calculations of the uni- 
versal critical exponents for this transition {i/ = 1.0 ±0.2 
and (3/v = 0.80 ± 0.05) show that it is indeed distinct 
from the Heisenberg transition. Direct simulation of the 
NCCP"'^ model however yields exponents that are consis- 
tent with these. Moreover, the large « 0.6 of the vector 
(magnon) field implied by these exponents can be heuris- 
tically understood since the magnon can decay into a pair 
of unconfined spinous at this critical point. 

Thus, an important conceptual issue that is clarified 
by our work is the existence of two different spin rota- 
tion symmetric critical points in D=2-|-l. The first is 
the Heisenberg transition whose 'soft spin' field theory 
is given by the 0(3) (f)^ theory. This describes of course 
the transition in the 0(3) sigma model if it is regularized 
by putting the spins on the lattice. It also describes the 
transition in the lattice CP^ model with a compact gauge 
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field. The second transition is described by a 'soft spin' 
field theory with a pair of complex scalar fields (that 
transform as spinors under spin rotations), coupled to 
a noncompact gauge field. This describes the ordering 
transition of the lattice 0(3) sigma model with hedgehog 
suppression. It also describes the transition in the lattice 
CP^coupled to a noncompact gauge field. Earlier indis- 
criminate use of what is effectively the NCCP^model to 
describe the Heisenberg transition are therefore to be re- 
considered. 

We also studied various physical extensions of the 
hedgehog suppressed 0(3) model, such as adding Zee- 
man field and also the effect of finite temperature. The 
P* phase survives in small external magnetic field and 
retains the photon, which is an evidence for gapped 
spinous. At finite temperatures, a thermal Coulomb 
phase with long ranged power law electric field corre- 
lations exists, and undergoes a Kosterlitz-Thouless tran- 
sition to the usual paramagnetic phase at higher temper- 
atures; this is indirect evidence for gauge charge carrying 
particles (spinous) that interact via the Coulomb poten- 
tial which is logarithmic in 2D. 

An interesting extension is obtained by breaking the 
full rotation symmetry down to the easy-plane XY sym- 
metry. We showed that this model possesses a remarkable 
'self-duality' property. Indeed as far as we know this is 
the first purely bosonic model in D=2-|-l that displays 
this property. In particular the critical point describ- 
ing the ordering transition in this model is self-dual. An 
important role is played by this critical theory in a forth- 
coming publication 19J, where it is conjectured that this 
critical point may already have been seen in numerical 
experiments on easy plane spin half quantum antiferro- 
magnet on the square lattice where surprisingly a contin- 
uous transition is observed between a spin ordered state 
and a valence bond solid that breaks lattice symmetries 

In D=2+l the P* phase and its gapless photon excita- 
tion only exist in the limit of infinite hedgehog suppres- 
sion. Thus, within this phase, the complete suppression 



of free hedgehogs represents an extreme limit that may 
seem unnatural. However, reasoning in this limit can be 
conceptually powerful and throw light on many tricky is- 
sues. Besides, we note that even if hedgehogs eventually 
proliferate and lead to a gap for the photon, this P* phase 
may be relevant for describing physical crossovers in some 
strongly correlated systems at energy scales above the 
photon mass. In contrast to the P* phase which in the 
CP^ language is unstable to turning on a finite monopole 
fugacity, the stability of the critical points where gapless 
gauge charged particles are present is a more involved 
question, and depends on the number of gapless fields 
present and the particular monopole creation operator 
under consideration. In fact, it is argued in [l^ that 
quadrupled monopole operators (which are relevant to 
the spin half quantum antiferromagnets on the square 
lattice) are irrelevant at the NCCP"'^ critical point in both 
the isotropic and easy plane limits. This would allow for 
a continuous transitions between valence bond solid and 
Neel phases in such systems, and these transitions would 
be controlled by the critical points studied in this paper. 
Finally, we expect a similar P* phase in the (3-f 1)D 0(3) 
model, which will now be present also with strong but fi- 
nite monopole suppression, which would provide a spin 
rotation invariant (3+l)D model which exhibits a frac- 
tionalized Coulomb phase with deconfined spinous and a 
true photon excitation. 
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